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Electron-Optical Shadow Method of Magnetic-Field 

Mapping l 

L. Marton, J. Arol Simpson, and S. H. Lachenbruch 

The essential features of the electron-optical shadow method for the quantitative 
mapping of magnetic fields are presented in unified form. The processes of obtaining and 
analyzing the shadowgrams are given in some detail. Iterative procedures of analysis, 
which have proved successful in practice, are illustrated by the mapping of an otherwise 
inaccessible ferromagnetic domain field. Sections 5 and 6 give an analytical method, involving 
inversion of integral operators. The method, although it has not yet been applied to the 
analysis of shadowgrams, is given because of its generality and applicability to the equivalent 
problem of interpretation of fluid dynamic interferograms. 



I. Introduction 

A technique was introduced about 3 years ago 2f 3 
for quantitative mapping by means of an electron- 
optical system, of electric- and magnetic-field 
distributions, which, because of their small dimen- 
sions extending to below 1 micron or because they 
would appreciably be affected by the presence of a 
probe, are not capable of satisfactory measurement 
by conventional probe methods. Since that time, 
both experimental and theoretical studies of this 
"shadow method" have been made. The purpose 
of this paper is to summarize briefly in unified form 
the essential features of the method, as modified and 
broadened through subsequent published and un- 
published studies and experiments. 

The shadow method may be considered as con- 
sisting of the following three main basic steps : [to be 
discussed separately] : 

Step 1 (Experimental). The field to be studied is 
placed in the path of an electron beam in the object 
space of an electron-optical lens system whose 



i This paper is a summary of 3 years work in the Electron Physics Section of 
the National Bureau of Standards. Many members, past and present, have 
contributed to the work described. The bulk of the mathematical development, 
part of which appears as sections 5 and 6, is due to S. U. Lachenbruch. In the 
experimental part, much is due to the work of A. Van Bronkhorst, G. Ford, 
J. Suddeth, and T. McCraw. 

2 L. Marton and S. H. Lachenbruch, J. Research NBS 43 409 (1949) RP2033. 

3 L. Marton and S. H. Lachenbruch, J. Appl. Phys. 20, 1171 (1949.) 



image space contains an opaque object, such as 
a wire grid, and the distortion of the shadow of this 
object due to the presence of the unknown field is 
measured as a function of position. When the field 
does not have axial symmetry, additional experi- 
mental data of value may be obtained by repeating 
the experiment with different orientations of the 
field and /or varying beam energies. Figure 1 shows 
the arrangement of one such nonsymmetrical sample 
and a representative shadow graph. 

Step 2 (Geometrical) . From this shadow distortion, 
and from known geometrical parameters of the 
electron-optical system, the angular deflections of 
different rays by the given field are calculated. 

Step 3 (Analytical). The field distribution is 
calculated from this deflection function and from 
any known or assumed qualitative properties of the 
unknown field. 

2. Step 1 (Experimental) 

Figure 2 is a schematic diagram of the electron- 
optical arrangement used. The main components 
of the system are : 

(a) an objective-lens system, characterized by the 
principal planes, P, P' , and the focal length, /; 

(b) the unknown field, placed in object space with 
its center at a distance s', from the first principal 




Figure 1. 

(a) Typical arrangement for obtaining the shadowgraphs of a nonsymmetrical object, (b) Typical distorted mesh shadows. 
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Figure 2. Schematic drawings of the optics of the shadow method. 



plane (this determines the conjugate distance s) ; 

(c) a point source, E', of electrons, at an "object 
distance" e' ', (usually e'=<x>) } which determines a 
crossover point, E, at a distance e (usually e=f) in 
the absence of the unknown field; 

(d) a wire mesh or other shadow-casting obstruc- 
tion, placed in image space at a distance c, thereby 
determining the conjugate distance c' ; 

(e) an optional projector lens, M, for magnifying 
the shadow image, characterized by the lens separa- 
tion, s + l, and the focal length, /i; 

(f) a fluorescent screen or photographic plate in 
the final image plane (S or Si) conjugate to the plane 
containing the center of the field. 

Most of the work at the Bureau has been done with 
a modified 4 electron microscope, which provides the 
above components in a convenient form. 

As the following conditions are to be assumed in 
the calculations, it is necessary that they be satisfied 
experimentally : 

Assumption 1. The trajectories within the field 
lie entirety in one plane (the plane of fig. 2) . 

Assumption 2. All lateral distances and angles, 
measured from the optical axis, are sufficiently small 
to make geometrical aberrations negligible, so the 
Gaussian dioptrics is applicable. 

Assumption 3. The distance between deflecting 
field and lens field is sufficiently great so that the 
interaction between the two fields may be neglected. 

Assumption 4- The electron energy, V, is in the 
nonrelativistic range. 

The assumptions place limitations on the form, 
position, orientation, and order of magnitude of the 
electromagnetic field. Thus, assumption 1 requires 
that planes exist in which the field produces only 
plane deflections, and assumption 2 requires that 
these deflections be small. No restriction is placed 



« J. A. Simpson and A. Van Bronkhurst, Rev. Sci. Instr. 21, 669 (1950). 



on the time variation of the field which may be 
static or varying. 5 

The solid ray in figure 2 represents a typical 
electron trajectory, deflected by an angle <t> by the 
unknown field before entering the lens system, and 
casting, on the final screen Si, a shadow Ai of the 
point G of the opaque obstruction. The broken 
ray represents the trajectory, which, in the absence 
of the deflecting field, would pass through this same 
point G, the corresponding shadow being at A i. 
The distance between these two points, which 
represents the displacement of the shadow of G due to 
the presence of the field, may be measured for each 
point G, whose shadow can be identified and then be 
plotted as a function of the distance a(SiAi), which 
is appr oximately proportional to the coordinate 
a'(S'A') in the deflecting field. 

Practically, the distances Aa can be measured 

(1) directly on a plate on which two exposures with 
and without deflecting field are superimposed; 

(2) if only one exposure with deflecting field is taken 
and there exists an undistorted region of the shadow 
pattern, the straight line corresponding to this part 
of the plot (a— do) versus a can be extrapolated, and 
the distances between the actually observed portions 
of the shadow and their extrapolated line can be 
measured; (3) if no undistorted part of the shadow 
pattern exists, a second mesh can be introduced 
between the electron gun (or condensor lens) and 
the deflecting field providing a reference mesh-wire 
system. 6 

Tests with the known field in the equatorial plane 
of a Helmholtz coil as deflecting object have shown 
that the errors introduced in this (a—a^) versus a 
plot by the aberration of the lens and the experi- 
mental uncertainties are smaller than 2 percent. 



5 L. Marton and D. L. Reverdin, J. Appl. Phys. 6, 617 (1950). 

6 D. L. Reverdin, J. Appl. Phys. 22, 257 (1951). 
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3. Step 2 (Geometrical) 

Assumption 2 above implies the existence of a 
linear relation between Aa h a measurement in the 
final image plane, and the angle <£, through which the 
unknown field deflects the ray that passes through 
the typical point A' of the field. This relation may 
be expressed as 



<f)=KAai 



(1) 



where K is a known function of the magnifications 
and other fixed geometrical parameters of the system. 
From the geometry of figure 2, it is easy to show 
that K may be expressed as 



K- 



1 



/ Ml (,--l)(l_ e / M ) 



(2) 



where 



H 1 =l 1 /l=AailAa= magnification of projector system; 

jj L =s/s f = Aa/Aa / = magnificat ion of magnetic object 
on screen S; 



v=(s — e)/(c—e)-- 



- magnification of opaque object C 
as determined by its shadow on 
screen S in absence of deflecting 
field ; 

~e/e f which vanishes when incident 
beam is parallel. 



Whereas in step 1 the deflecting effect of the field 
was magnified, by means of a lens system and 
shadow-casting object, in order to facilitate the 
measurement of that effect, in step 2 this magnified 
effect, represented by the measurement Aa iy is con- 
verted back to angular deflection in object space. 

The problem of experimentally determining K by 
measuring the quantities on the right-hand side of 
eq (2) is somewhat complicated by the distortion 
introduced by the projector lens. This distortion 
makes the precise determination of the magnifica- 
tions difficult. A series of tests 7 on the known field 
of a small Helmholtz coil shows that with care the 
errors introduced in the determination of K can be 
less than 5 percent. 

4. Step 3 (Analytical) 

The discussion of this step will be limited, for 
simplicity, to magnetic fields, which have been 
studied in considerable detail. 

We note first that assumption 1 above confines the 
analysis of magnetic fields to planes (such as the 
equatorial plane) that intersect all magnetic lines of 
force orthogonally. All electrons initially moving in 
such a plane will remain therein because their deflec- 
tion by the field is everywhere normal to both the 
path and the lines of force. The plane of figure 2 
must therefore have this property. 

7 J. A. Simpson, Proceedings of the NBS Semicentennial Symposium on 
Electron Physics, NBS Circular 527 (1954). 



Figure 3, which lies in one such plane, summarizes 
the theoretical problem of step 3. The polar coordi- 
nate system (r, 0) is attached to the magnetic field 
H(r, 6), which is directed normal to the plane of the 
figure. The relative orientation of field and optical 
system may be varied experimentally. The curve L 
represents a typical electron trajectory through the 
field in this plane. It may be considered as approxi- 
mately straight because assumption 2 implies that 
the deflection angle <j> is necessarily small — otherwise 
the method, which involved geometrical optics, would 
not have been applicable. Moreover, this assump- 
tion eliminates the variation of beam energy as an 
added parameter for investigation of nonsymmetrical 
fields. 

In steps 1 and 2 a lens system has served the pur- 
pose of measuring, by magnifying its effect, the 
deflection angle, <£, corresponding to each of a certain 
set of directed lines, L, in the plane of figure 3. 
Now, each such line is identified uniquely by a pair 
of coordinates (j) y f$) as shown, where the angle p 
may be varied experimentally by arbitrary plane 
rotations of the field relative to the lens system, and 
p is the classical impact parameter, and in the case 
of parallel illumination is equal to a /ju. Having 
served its purposes, therefore, the lens system may 
now be disregarded, and the problem remaining 
may be stated as follows: For each of a certain set of 
directed lines L=L(p,(3) in the plane of figure 2, each 
identified by coordinates (p,P), we have an experi- 
mentally determined deflection angle <j>=<j>(<pfi) y and 
hence, by the well-known relation (giving the 
deflection of a charged particle in motion through 
a magnetic field), the value of 






(3) 



the integral of H along that line. (Here e, m, and 
V are the charge, mass, and potential of the elec- 
trons.) We wish to derive from this the magnitude 




Qirfi) 



Figure 3. Coordinates used in the 
analysis of shadowgraphs. 
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and sign of the unknown H at an arbitrary point Q 
with polar coordinates (r,d). It is in this step — 
that of obtaining the unknown function H (r,0) 
from its integral I{pfi) — that a number of procedures 
have been proposed, and several tried. These 
methods are summarized below. 

It is worth noting that the problem is essentially 
the same one that arises in the interpretation of 
interferograms taken of air flow, and any of the 
methods used in that field apply. Recently, it has 
been shown that the problem may be programmed 
for a large-scale computer. 8 If such a computer is 
available, this represents the ideal method of 
analysis of shadowgrams. 

For the case where the field is known to have 
axial symmetry about some axis normal to the 
plane of figure 3, we can choose the origin as the 
intersection of the axis with the plane. Then 
H=H(r) is independent of 0, and also I=I(p) is inde- 
pendent of j8, so the problem is greatly simplified. 
Some of the methods about to be described are 
limited to this case, which, due to the difficulties of 
the more general case, is the only one that has been 
investigated in detail up until this time. 

4.1. Method (a). Fitting of parameters 

In the method first used, the unknown field 
H{rfi) is assumed to have a specific analytical form, 
usually selected from known or predicted physical 
properties of the field and field-producing object, 
and involving one or more unknown parameters 
{ai,a 2 , . . .) to be determined 



H(r,0)=F(r,6,a l9 a2 . . . a n \ 



(4) 



where F is a preassigned function of r and 0, periodic 
in 0. 

Integrating the specified function F in advance 
gives a predicted analytical form for /, as a known 
function G of (n+2) quantities 

I(P,P)=\ F(r,6,a 1 ,a 2 . . . a n )ds = G(p ) p j a 1 . . . a n ) 

(5) 

We wish to solve for the n quantities a t by matching 
this (7-f unction against the values of the /-function 
obtained experimentally for different values of p 
and j8. 

One way to do this is statistical- — by applying least 
squares to the experimental results, starting with an 
equation of the form 



TO 



■G(p k ,p k ,a,i . . . a n )] 2 = minimum, (6) 



m being the number of electron paths for which 

8 F. D. Bennett et al. Interferometric analysis of airflow about projectiles in 
free flight, Ballistic Research Laboratories, Aberdeen Proving Ground, Report 
No. 797 (March 1952). 



experimental data have been obtained. Here the 
closeness of fit can usually be used as a partial check 
on the validity of the functional form F chosen. But 
how good a check is not known until we consider the 
question of uniqueness, mentioned at the conclusion 
of section 6. 

An example of the method of fitting of parameters 
is described in detail in an earlier paper (see foot- 
note 3). The field of a magnetized wire, in its 
equatorial plane, was evaluated on the hypothesis 
that its form approximated that of the fielol between 
two equal and opposite point poles of unknown 
magnetic moment, a 1 , separated by an unknown 
distance, 2a 2 . For this case, the function F of eq 
(4) is 

H(r) = F(r,a l ,a 2 ) = a 1 (r 2 +ai)-^. (8) 

The corresponding function G is found from eq (5) 
to be 

G(p, a lf a 2 )=2a 1 /(p 2 +a%). 



(9) 



The moment a x and pole separation, 2a 2 , and hence 
the field H(r), were obtained by a least-squares 
method of fitting the function G to an experimental 
scatter diagram derived from deflection measure- 
ments. 

4.2. Method (b). Evaluation by step-function 
approximation 

In a second method, the field is assumed to exist 
only in a finite portion of the plane of figure 3, and 
this portion subdivided into n regions B k of chosen 
shape and size, and the field calculated as if it were a 
step-function constant within each such region 



H(;r,d) = a k 



(10) 



where (r,0) eE k ,k=l,2 . . . n. By integration we 
obtain 



fc=i 



(id 



where \ k (p,P) is defined to be the length of the 
segment of intersection of the line L(p,fi) and the 
region R k . 

If, then, / (p,P) is known experimentally for n 
different electron paths L(p,p), there result n 
simultaneous linear equations in the n constant 
values a k of H. These equations may, in general, be 
solved for the a k , although the labor may be great in 
many cases. Moreover, it is not possible generally 
to choose a grid and a set of rays that insure that the 
set of equations is linearly independent. This 
proves to be a very serious difficulty. 

The technique to be described avoids both diffi- 
culties. Figure 4 shows the grid used to study the 
field about a corner of the cobalt crystal of figure 1. 
This is merely an illustration, and a number of other 
grids have been developed for specific problems, 
but all use the same technique. The rays are so 
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HEXAGONAL FACE 

COBALT CRYSTAL 



Figure 4. Analysis of magnetic fringe field of a cobalt crystal. 

(a) Grid used in evaluation of fields about the comer of the sample by means of the step-function approximal ion . (b) Results of this method on the fringe field 
of a single cobalt crystal. 



chosen that each ray passes through only one square 
for which a k has not previously been determined. 
The first three such rays are shown in figure 4. With 
such a system, the determinant is triangular and can 
be solved within a reasonable time. Unfortunately, 
with a given number of values of the angular param- 
eter, only a limited number of grids may be used 
before it is impossible to draw a line that does not 
pass through two or more unknown grids. In the 
grid illustrated, this takes place after 47 squares 
have been evaluated, when there are shadowgraphs 
available for each 15° of revolution. Figure 4, b, 
shows the results of this method applied to the field 
of the cobalt crystal 

The errors inherent in this procedure are compli- 
cated and large. In addition to any experimental 
errors involved, one has the error caused by replacing 
the integral with the rather crude sum. This error 
can be estimated by changing the number of the 
meshes up to the maximum possible and observing 
the change in values. A value for this error seems, 
in practice, to be as large as several hundred percent. 
The exact value depends on the magnitude of the 
gradient and also on the sensitivity of the original 
shadowgraph. (The field is considered to fall to zero 
outside the grid, which, in practice, is extended to 
the limit of the region in which there is a measurable 
deflection.) 

Known properties of the field sometimes suggest a 
special choice of subdivision. For example, if the 
field is known to have rotational symmetry about 
some known point in the plane of figure 3, the most 
appropriate mode of subdivision would be into rings 
bounded by concentric circles about that point, 
inasmuch as the field is constant along any one of 
these circles. 



The field surrounding a magnetized wire has been 
computed by this method and showed good agree- 
ment with the result obtained by method (a), in 
which the field was assumed to be that of a pair of 
point poles (see eq 8 and 9). 

4.3. Method (c). Series expansion 

A field that is known to have axial symmetry 
about a point at which it may reasonably be assumed 
to become infinite can, in general, be assumed ex- 
pansible as a series in inverse powers of the distance 
r from that point, such a series being convergent for 
sufficiently large r 



a 2 . a A 



ff('H3+3+ 



(12) 



Physically, this expansion corresponds to assuming 
that the field is that of an //-pole, and hence no terms 
below that corresponding to a single pole (a z ) appears. 
By integration of (12) along the horizontal line 
L(p, t/2) whose equation is r sin 6=p, we obtain 

I(p)= f H(r)ds= f H{r)d{J^tf) 

JL{p,irl2) JL 

=2/7^+...)^= 



_ira 2 , 2a 3 



(13) 



On matching this series expression for I(p) against 
the function obtained experimentally, one is able to 
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evaluate numerically the coefficients a 2 , « 3 , . . ., and 
hence deduce the series representation (12) for H. 

Even when the condition of axial symmetry is not 
satisfied, it may sometimes be possible to find types 
of series appropriate for expansion of given field 
distributions. The numerical evaluation of the co- 
efficients would then be obtained, as above, by 
matching the termwise integral of the series with the 
experimentally derived function I(p,P) . The restric- 
tions, techniques, etc., would of course vary from 
case to case; thus, for example, if a Fourier series, or 
other series in orthogonal functions, is used, the 
unknown coefficients will generally be obtainable 
from an explicit formula involving definite integrals. 

The series expansion method, although it has a 
certain amount of generality, must be used with 
some caution. On account of questions of con- 
vergence, for example, a series of form (12) could 
never be used in the neighborhood of the point 
r=0. 

5. Direct Method of Inversion of Integral 
Operators 

It is to be noted that in all three methods described 
above, the approach consists in assuming a priori 
some specific analytical form for H, involving one 
or more unknown parameters, which are then eval- 
uated from experimental results. In this sense, they 
may all, in fact, be regarded as special cases of 
method (a). If, for example, the function F of (4) 
is taken to be the step function of (10), the result 
is essentially method (b) ; and if F is a series such 
as (12), one obtains method (c). Still other methods 
ma} 7 be derived by choosing F as a function of still 
different form. Thus a generalization of method (b) 
is obtained by choosing F to be a linear or second- 
degree function, rather than a constant, in each 
region R k . 

An entirely different approach consists in inter- 
preting (3) as an integral equation, and observing 
that it can, in general, be uniquely solved for the 
unknown function H(rfi) in terms of the function 
7(p,/3) which may be obtained experimentally. 

To date, no use has been made of the method in 
experimental work. However, because of its gen- 
erality, preciseness, and application to the analysis 
of interferograms in fluid dynamics, it is given in 
detail. 

5.1 Analytical solution for case of axial symmetry 

For the following considerations, the integral 
equation (3) may be expressed in a more specific 
form by using the equation of L(p,I) (see fig. 3). 

P=r cos (0-0). (14) 

By solving (14) for 0, we see that, on the line L(p,P) 

(/3 + cos -1 p/r on the one half-line (15) 



0= 



j8 — cos -1 p/r on the other half line 



Substituting (15) into (3), and expressing 
function of r, we obtain 



as a 



(16) 



We take first the important special case where 
the field is known to have axial symmetry about 
some axis normal to the plane P. We can choose 
as origin of coordinates the intersection of this axis 
with P. Then H=H(r) is independent of 0, and 
also I=I(p) is independent of 0. Hence, (16) sim- 
plifies to the form 

I(p)= 2 I " Htr)rdrU!?^p~ 2 . (1 7) 

Jp 

Now Abel's integral equation, 

m= V g(y)dyl(z-yY (0<x<i), (18) 

J a 

is known to have the unique continuous solution 



g(y)= 



sin Xtt d C y 



dy 



f(z)dzl(y-zy- 



(19) 



provided f(z) satisfies certain very general continuity 
conditions and /(a) = 0. If in (18) one replaces z 
by (—P 2 ), /(z) by I(p), a by - oo , X by \ by 2, y by 
(— r 2 ), and g(y) by H(r), it is seen to reduce exactly 
to (17). By (19), therefore, the unique solution of 
integral eq (17) is 



H(r)- 



sin tt/2 d 



d(-r 2 ) 



§~_j(v)d{-v 2 )Nv 2 



which may be reduced to any of the following- 
equivalent forms 

irr Jo dr 

=_I fVcvF+^MAFH? 

TTJO 

We have thus an explicit expression for the field 
H(r) in terms of the empirically known function I{p). 
If I(p) is known at a sufficient number of experi- 
mental points so that a curve may be plotted, the 
value of H(r) corresponding to each value of r may 
be calculated directly from (20) by some numerical 
integration procedure. 
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The integral eq (17) and its solution (20) may bo 
abbreviated as 



I(p)=T{H(r)}; H(r) = T^{I(p)} 



(21) 



This form suggests a somewhat different way of 
dewing the problem: viz, as the problem of inverting 
a linear integral operator T in function-space. 

It can also be considered as a problem in integral 
transforms. When miscellaneous elementary func- 
tions are substituted for H(r), and the corresponding 
functions J(p) are calculated from (17) and simpli- 
fied, and then verified by substitution into (20), a 
table of transforms may be constructed, similar in 
form to a table of Fourier or Laplace transforms, 
and the linearity property may be used to derive the 
transforms of more complicated functions, including 
the field-distribution functions most frequently 
encountered. 

The one-to-one correspondence implied by (21), 
or by (17) and (20), provides a conclusive answer, 
for the axially symmetric case, to the question of 
uniqueness of the result — a question not rigorously 
investigated previously. We can now state that if 
the shadow method is applied to a magnetic field 
known to have axial symmetry, aid if a sufficient 
number of experimental measurements are made so 
that the empirical function /(/;) may be considered 
as a known continuous function of the continuous 
variable p, then there 4 exists one and only one 
(continuous) function H(r) satisfying (17), and it is 
given in terms cf the known I(p) by (20). 

5.2. Analytical solution for fields of general form 

Finally, we turn to the general case, removing the 
condition of axial symmetry. The general problem 
is that of solving (16) for the function H(rfi) of two 
variables, when sufficient experimental data are 
available so that I(p, f$) may be treated as an experi- 
mentally known function of both variables. 

We note first that for any fixed value of r, the 
field H must be a periodic function of 6, with period 
2tt, and so it must have a unique Fourier series 
expansion with functions of r as coefficients 



H(r,6)= S H n (r)exp(inO). 



(22) 



Similarly, / has a unique Fourier expansion, with 
functions of p as coefficients 



^(p,i8) = Z)A(p)exp(iwj8). 



(23) 



On substituting (22) and (23) into (16), and equat- 
ing coefficients of like powers of exp (if$) on the two 
sides of the resulting equation, one obtains 



J p 
(n - 



I„(p) = 2\ H n (r) cos (n cos l p/r) Z-^-= > 



(24) 



For each value of n, (24) is an integral equation for 
H n (r) in terms of I n (p). 

A calculation carried out in the appendix shows 
that these integral equations have unique solutions 
for the //-coefficients in terms of the corresponding 
/-coefficients. If, then, one has sufficient experi- 
mental data to define the continuous function I(p,f$), 
the coefficients in its Fourier expansion (23) are 
uniquely determined from the formula 



In(p> 



Z7T Jo 



(25) 



the coefficients H n (r) are then determined in turn as 
solutions of the integral eq (24), and so the desired 
function H(rfl) will be obtained explicitly in the 
form of a Fourier series (22). 

The fact that a Fourier series representation was 
used here as a tool to solve the general problem, and 
to investigate existence and uniqueness questions, 
may tend to obscure the basic nature of the result. 
By analogy with (21), we may abbreviate the integral 
eq (16) and its solution as 

I(]>fi) = T{H{r,9)\ and H(r,9) = T- l (I(pfi)). (26) 

This form again suggests that the problem may be 
viewed as that of inverting a linear integral operator 
T, this time in a space of functions of two variables. 
It is possible that more compact expressions and 
proofs might be possible when the problem is attacked 
from such a point of view. 

The question of uniqueness is again seen to be 
conclusively answered as a by-product of this anal- 
ysis — this time for the general (nonsym metric) case. 
We can state that if, in applying the shadow method 
to a magnetic field that is not known to have axial 
symmetry, sufficient experimental data are obtained 
so that the function I(p,P) may be considered as 
a known continuous function of the two continuous 
variables p and |3, periodic in f$ (with period 27r), 
then there exists one and only one (continuous) 
function Z/(/',0), periodic in 6 (with period 27r), and 
satisfying (3). It is obvious that in order to satisfy 
the requirement of sufficient data in this general 
case, it will be necessary to determine the dependence 
of /on p by repeating the experimental measurements 
for different relative orientations of field and optical 
axis. 

6. Appendix. Outline of a Procedure for 
the successive solution of the integral 
equations (24) 

(a) When n = 0, (24) becomes 

/ (p)=2f 00 i/o(r)dr/Vr2^o (27) 

Jp 

which has exactly the same form as (17), and, therefore, by 
(20), has the unique continuous solution 



//o(r) = --P/o(p)/Vp 2 - 

Wjr 



(28) 
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(b) When n=±l, (24) takes the form 

Jp 
or 

Jp r ^ r 2 — p 2 



(29) 



which is of the same form as (27) except that the functions 
h(p) and H (r) are replaced, respectively, by I±:(p)lp and 
H±i(r)/r. The unique solution, obtained from (20), is 
therefore 



irjr dpi p J J V 2_ 



Vp 2 — r 2 



(30) 



(c) The higher-order coefficients (for |n|>2) are more diffi- 
cult to evaluate. The following method of obtaining the 
coefficients for n = ± 2 and n = ± 3 appears capable of extension 
to coefficients of higher order. 

Since 

cos[±2 cos -1 k] 

= 2(/c 2 -l)-f 1; cosIiScos" 1 k] = 4k(k 2 - 1) + jfc (31) 

(24) becomes, for n= ±2 or ±3, 

where 

A 2 = A- 2 =2, A 3 = A- 3 = 4 (33) 

But (32) may be written as 

~~ T^n:o = ^ , i o , =^— 2A n — r ^- yr 2 — p 2 dr, (34) 

plnl-2 J p rl n| - 2 y r 2_p2 Jp f\n\-l * f ■ 1 \ J 

which becomes, on applying integration by parts to the last 
integral, 

where 



y r 2 _ ^2 



f (v ^ H n(r) A C^Hni^dp 



(36) 



r |n|-2 Jr M , 

Equation (35) is again of the general form (17) with certain 



replacemen s, and hence, by (20) 

/n(r)= _ir ro i_r / 



i r d ri n (p)l 

-irjr dplp^\- 2 J 



dp 



Itfl- 



SV 



(37) 



Substituting (36) into (37) and differentiating with respect 
to r, 



py,(rn _,A n r H n (r) l -1 d r- d rz.^n 

|_ r ini-2j"t" r L r ui- 2 J 7T dr J dplp\ nl ~ 2 J 



dp 



-y]p 2 — 1 



(38) 



But (38) is a linear first-order differential equation in the 
unknown function i/ n (r)/r ln l~ 2 , which can be solved explicitly 
for H n (r) by the usual formula, and the arbitrary constant 
evaluated from the condition that it vanishes at r= <» . 
This determines the functions H 2 , H- 2} H 3 , H- z uniquely, as 
was to be proved. 

7. Summary 

The electron-optical shadow method of field 
mapping has been used at the Bureau for 3 years to 
study such diverse and otherwise inaccessible fields 
as the fringe fields of ferromagnetic domains of single 
crystals and collodial particles, space-charge fields 
in magnetrons, and the external fields of magnetic 
recording wire. The over-all accuracy of the method 
depends upon the analysis of the experimental data, 
and, in favorable cases, has been as good as 5 percent 
although 10 or 100 times greater error may result in 
the more unfavorable cases. But in any case, the 
method opens a new magnitude of fields to study 
and measure. 
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